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Topological way of thinking now goes beyond conventional solid materials, and topological char¬ 
acterization of classical mechanical systems governed by Newton’s equation of motion begins to 
attract much attention. To have a deeper insight on physical meaning of topological numbers in 
mechanical systems, we demonstrate the use of the Hannay angle, a classical counterpart of the 
Berry phase, as a symmetry protected topological order parameter. We first derive the Hannay 
angle using a canonical transformation that maps the Newton’s equation to the Schrodinger type 
equation. The Hannay angle is then used to characterize a simple spring-mass model topologically 
with a particular focus on the bulk-edge correspondence and new aspects of the symmetry in a 
classical system. 


A study of topological phenomena in solids has its ori¬ 
gin in quantum Hall systems |Tj. However, topological 
phenomena are not limited to quantum systems, but also 
found in classical systems such as photonic crystals gov¬ 
erned by the Maxwell equation Hi- Very recently, clas¬ 
sical mechanical systems obeying the Newton’s equation 
of motion begin to get much attention as a playground 
of topology SHU]- Although topology is an abstract 
mathematical concept, physical observables are charac¬ 
terized by the topological quantities and well described. 
For instance, in the quantum Hall systems, a topological 
number, i.e., the Chern number defined using the Bloch 
wave functions P2 , has a direct relation to the Hall con¬ 
ductance and the number of edge states m Another 
example is the Berry phase, which describes the topolog¬ 
ical edge states of Dirac fermions in solid state materials 
such as graphene [T4ffT7] . In contrast to the Chern num¬ 
ber, which is quantized by definition, the Berry phase 
is quantized and topological only with a help of some 
appropriate symmetry. It is the symmetry protection 
of topological phases. In any cases, nontrivial bulk as 
a topological phase implies existence of localized edge 
states if the system has boundaries, i.e., there is the bulk- 
edge correspondence. This is surely quite important for 
topological description of classical systems, where we do 
not mostly have bulk topological observables such as the 
Hall conductance. Even though the classical system is 
topologically nontrivial, we may only access its toplogi- 
cal character through the edge states. Direct access to 
the bulk topological observable in classical systems is an 
important issue to be addressed. 

As a first step toward better understanding of topolog¬ 
ical numbers in classical systems, here we demonstrate 
the use of the Hannay angle [iTS], which is a classical 
counterpart of the Berry phase, as a symmetry protected 
topological order parameter. We first present a concise 
formulation of the geometrical angle, i.e., the Hannay 
angle, in mechanical system using a canonical transfor¬ 
mation, such that the Newton’s equation is mapped to 


the Schrodinger type equation. Then, the solid relation 
between the quantized Hannay angle by symmetry and 
the edge states is established by numerically analyzing 
a simple spring-mass model. Interestingly, the symme¬ 
try protected nature is observed as boundary condition 
dependence of the edge modes. 

Let us begin with describing a geometrical phase in a 
classical system using a simple example, 2D motion of a 
mass point (mass m = 1 is assumed for simplicity) in the 
potential V,j,(x,y) = a 2t xT^x/2 with 

r 0 = 1 + A(1 — r) + 77 COS</>)<x 3 + A 77 sin , (1) 

where <f> has some time dependence and Oi is the Pauli 
matrix. V^(a :,y) is periodic in <fi with period 2-7T. This 
problem is so simple, but captures essence of the ge¬ 
ometrical phase in classical mechanics. For the fixed 
t f >, this system has two normal modes with frequencies 

uj± = Now, assum¬ 

ing the time dependence of <f> as fit, we follow the time 
evolution starting from the initial conditions x = ( ( 0 , 1 ) 
and x = 0 , which excites only the mode with toff. 
We assume slow evolution of <j>, i.e., fl/a -C A, and 
then, the other mode with wj is expected to be un¬ 
excited during the time evolution. In order to factor 
out rapid oscillation, which corresponds to the dynamical 
contribution, new dynamical variables x c = \x c ,y c ) and 
x s = \x s ,y s ) are introduced as x = x c cosd t + x s sm0 t 
with 0 t = fgCOQ t ,dt', and the equation of motion is di¬ 
vided into the coefficients of cos(9 t and sin# t . 

The numerically obtained time evolution of x c and x s 
is plotted in Fig. [T] We find that at t = T, when the 
potential gets back to the original form, x c {T) is same 
as ® c (0) for 77 < 0.5 while x c (T) is minus of x c {Q) for 
?7 > 0.5. Since the rapid oscillation is already factored 
out, this means that the system acquires extra 7 r phase 
other than the dynamical contribution for 77 > 0.5. The 
physical origin of this extra 7 r phase is actually simple. 
Namely, V^(x, y) is constructed by superposing two func- 
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FIG. 1. Upper panel: Contour plots of the potential V,f,(x,y) 
for y = 0.4 and 0.6. Lower panels: Time evolution of x c and 
x a with a = 2ir, A = 0.5, and T = 1200 for selected y. 


tions Vi(x, y) = cnx 2 +ci 2 y 2 and v 2 (x,y) = c 2 ix 2 +c 22 y 2 
with appropriate c. tJ , and the evolution from <f> = 0 to 
(f> = 2ir is achieved by relatively rotating v 2 (x, y) against 
v\(x,y) by 180°, instead of 360°. When v 2 (x,y) is domi¬ 
nant, 180° rotation results in a flip of x c after completion 
of one period, while when i’i(x, y) is dominant, such a flip 
does not take place. 

In order to interpret the extra n phase as a geomet¬ 
rical phase, i.e., the Hannay angle, we introduce action- 
angle variables [191, 20]. We start with a Hamiltonian 
H = (p' ■ p' + ^'T\q')/2 1 where p' and q' are canoni¬ 
cal variables, and Ta is a symmetric matrix describing 
the potential that is parameterized by A. Using the idea 
of the Bargmann variables |21| . we introduce a canoni¬ 
cal transformation by a generation function E(q',q) = 
— i( f q' — a^q^xiq 1 — a_q)/2 with a ± = \/2 ± 1 and 
7 a = Vt\, which leads to a transformed Hamiltonian 

H = i‘p 7 A q. (2) 

The “offdiagonal” form with y/l\ reminds us the “super- 
symmetric” treatment of the spring mass model proposed 
by Kane and Lubensky [B|. The equation of motion is 
now written as 

q = dH/dp = ij x q, P = -dH/dq = -i^fxP (3) 

with general solutions q(f) = e lA/xt qo andp(f) = poe~ lA/xt . 
Since Eq. ([ 3 ]) is exactly like the Schrodinger equation, 
i.e., it only contains the first order derivative in time, it 
is possible to give a “quantum” derivation of the Hannay 
angle. Namely, we can regard 7^9 and ip as a ket |^>) 
and a bra (ip\, respectively, and follow the time evolu¬ 
tion in a quantum fashion. (The variables q and p are 


independent, but in order to keep q' and p' real, they 
should fulfill the relation (ip)* = 7 a<Z-) Here, in order 
to emphasize the classical nature of the phenomenon, we 
derive it in a purely classical fashion. 

The action is defined as 


S{q,qo,Po, A) = / p • dq, (4) 

J C(q,q a ,paA) 

where C(q , qo,po, A) is a path in the phase space whose 
initial and final values of q(t) are q and qo, and initial 
value of p(t ) is p 0 , respectively. Explicitly, we have 


5(q,qo,p Q , A) = / p • qdt = ^ \og(Q a /Qo a )PoaQoa 

^ a 

( 5 ) 

with Q = Oq and P = pO' where O is a matrix diag¬ 
onalizing 7 a- Namely, 77 = O^uiO where 6j is diagonal 
with u> a as diagonal elements. Now, the action variables 
I a are introduced as I a = iPoaQoa- By definition, 2irl a 
represents increment of S against completion of one pe¬ 
riod of the normal mode with uj a . Using /, we can express 
S' as a function of q, qo, -f, and A. Then, we consider a 
canonical transformation taking S as a generation func¬ 
tion. Since I a is an adiabatic invariant, I a is an ideal 
quantity as new momentum. The new coordinate conju¬ 
gate to I a is 

<9S(q,q 0 ,p 0 ,A) . n (r , s 

= -WV- = -l\og{Q a /Q 0a ) = W Q t, (6) 


which is named as the angle variable. From this, the 
motion of mass point is obviously periodic in w a with 
period 27 t. Therefore, difference in ti 7 a indicates differ¬ 
ence in phase of oscillating motion. 

Next, we consider the case that A weakly depends on 
time. Then, the generation function S depends on time 
through the time dependence of A, and the Hamiltonian 
is transformed as H ^ H + (dS/d A)A. Using the trans¬ 
formed Hamiltonian, the equation of motion for ui a be¬ 
comes vj a = dH/dI a + d 2 S(q,I,X)/(dI a dX)X and we 
formally obtain 




dH 

dP 


d t H - 


&S(q,I, A) dA 
dl a dX d t 


( 7 ) 


The first term represents the dynamical contribution, 
whereas the second term represents the geometrical con¬ 
tribution, which is the Hannay angle 9 ^. Following the 
Berry’s argument [ 22 ] , 9 ^ is evaluated as 

e ’ = -i< 8 > 


for the case that A gets back to the original value after 
evolution. In practice, we should average out rapid oscil¬ 
lations to have a convenient formula. The best method 
for the averaging is applying integration over w a in the 
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range [0, 27 t] [T5]. Here, we focus on the situation that 
only one normal mode is excited. In that case, provided 
that a = a corresponds to the excited mode, I a is finite 
only for a = d, and it is reasonable to take integration 
only over We finally obtain 


n(H) _ d 

“"as 

d 




r 27T 


d tUo 
27r 




= ij(dd x d^) Sl5l &\. 


( 9 ) 


If we regard O a i as a wave function, Eq. § is exactly 
the Berry phase. 

Now, it is possible to give two interpretations for the 
flip of x at t = T for 77 > 0.5 in Fig. HI First, since 

is a symmetric matrix, Ox—^ in Eq. Q can be a real 
matrix, which results in 0^ H ' 1 = 0. However, in this case, 
it is impossible to make O$ continuous and periodic in <f> 
with period 2n. Namely, if we force O 0 to be continuous 
in cj), we have O 2 K = — Oq, which explains the flip of 
x. On the other hand, it is possible to amend so 
that it becomes continuous and periodic by introducing 
a complex phase factor, typically e 1 ^ 2 . Such a phase 
factor leads to C> 2 tt = Go, but then, Eq. § gives an extra 
phase factor of n, which explains the flip of x. Possibility 
of two kinds of interpretations corresponds to the gauge 
degrees of freedom in quantum cases. 

This argument indicates that 9^ should be 0 or 
7T, since = ±Oo is the only possibility with 

the gauge choice in which O$ becomes real. How¬ 
ever, 9 I ' H ' ) deviates from 0 or 7 r if the starting Hamil¬ 
tonian is changed to break the time reversal sym¬ 
metry. As an example, let us consider the situa¬ 
tion that the system is rotated with angular frequency 
Hr. Then, the Hamiltonian written in the rotating 
frame is obtained by replacing p’ by pi + if! RO y q[. 
Again, we follow the time evolution of x c and x s in¬ 
troduced as q[ = x c cos9 t + x s sin9 tl but with modified 

= a\J 1 + /3 2 ± \J A 2 (l - 277(1 - 77 ) ( 1 - cos^)) +4 (3 2 
where /3 = Hr/ a. We use the initial condition x = 
*(0,1) and x = *(uo,0) with vq = —2/3(1 + f3 2 — 
y/A 2 + 4/3 2 )/(A - 2/3 2 + y/A 2 + 4/3 2 ) in order to excite 
only the mode for to// . 

The numerically obtained (x c , x s ) for 77 = 0.4 and 
77 = 0.6 are plotted in Figs.[2](a) and[2jb). In contrast to 
Fig. [I] where x s and y s are negligible for all time, x s and 
y s have significant time dependence. Furthermore, we 
have x c (T ) 7 ^ ±cc c ( 0 ), which indicates the extra phase 
factor other than 0 or 7 r. With our choice of the initial 
condition, the extra phase factor 9 corresponding to 9^ H ' t 
is evaluated as 9 = arccos y c (T), and plotted in Fig. [ 2 ] as 
a function of /3 = ttii/a. We see that 9 clearly deviates 




FIG. 2. (a,b) Time evolution of x c and x a for 77 = 0.4 and 
0.6 for /3 = 0.1. (c) /3 dependence of 6 for y = 0.4 and 0.6. 
The parameters a = 2n, A = 0.5, and T = 1200 are used in 
(a)-(c). 


from 0 and n for finite Hr. Now, it is desirable to find out 
a transformation mapping the rotated Hamiltonian to the 
Hamiltonian like Eq. ([ 2 ]). However, a naive extension of 
our transformation by F(q',q) is only applicable in the 
case that [cr^F^] = 0, which limits us to the trivial case 
with r 0 oc 1. A possible analytical formulation for the 
rotated case is a future problem. 

For the demonstration of the usage of the Hannay an¬ 
gle as a symmetry protected topological order parameter, 
we introduce a one-dimensional “dimerized” spring-mass 
model schematically depicted in Fig. [3](a). The involved 
parameters are mass of the mass points m = 1, frequen¬ 
cies associated with sublattice dependent harmonic po¬ 
tential wo±, and spring constants k\ and k 2 . For later 
convenience, we introduce parameters u>q, uq, uq, rj as 
w o± = co>q ±ct;§, k\ = 77 w 2 and k 2 = (1—77)w 2 . We limit the 
motion of mass points along the ID chain, then, the dy¬ 
namical variable for this system is deviation of each mass 
point from the equilibrium position, which is written as 
x na for the mass point at nth unit cell and sublattice a. 
Assuming the periodic boundary condition and introduc¬ 
ing new variables fka as x na = e~ mk f ka , the equation 
to be solved becomes f k = - H k f k where f k = %f k i,f k2 ), 
and 


Hk = (wq + w i)i - w 2 Reg fc <7i - w 2 Im g k a 2 + u%a 3 (10) 


with g k = (1 — rj) + ye lk . Now, H k contains complex el¬ 
ements, but it is straightforward to extend the former 
arguments for the case that F is real symmetric }23j . 
For fixed k , there exist two eigenmodes with frequencies 

uj±(k) = + w 2 ± y/wj + w 2 |g(fc)| 2 . When oj 2 = 0, 

H k does not contain a 3 and then, we say that the system 
has a “chiral symmetry” since H k — (wq + w 2 )1 anticom¬ 
mutes with (J 3 . Although it is not the chiral symmetry in 
a strict sense because of the constant diagonal part, this 
symmetry is important in topological characterization. 


Furthermore, for = 0, it is possible to map Eq. (10) 
to Eq. (JT]) by a parameter independent transformation 
such that <7i —> &3 and 02 —> a 1 , which indicates that 
the two problems are essentially identical to each other. 
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FIG. 3. (a) Schematic picture of the “dimerized” spring-mass 
model. An employed unit cell convention is also shown, (b) 
0 J 2 dependence of 9 in the spring-mass model for 77 = 0.4 and 
0.6 with loq = wi = 27 t and T = 800. (c,d) Frequency spectra 
for the finite size spring-mass model with the (c) fixed and (d) 
free boundary condition. (e,f) The fixed and free boundary 
conditions. 


Next, we follow the time evolution regarding k as a 
controllable parameter driven as k = f H. More specifi¬ 
cally, we solve 

^ = -*»,/«). (in 

with the initial condition /( 0 ) = fo = %1/V2, l/\/ 2 ) and 
/( 0 ) = 0 , which excites only the lower branch with w_( 0 ) 
at t = 0. In order to factor out the rapid oscillation, i.e., 
the dynamical contribution to nj Q , new variables u±(t) 
are introduced as f(t ) = e 10t u + (t)+e~ l6t U-(t) with 9 t = 
/ oWQ r dT, and we solve equations for u±(t). The initial 
condition (/( 0 ), /( 0 )) = (/o,0) is satisfied by the choice 
of u±( 0 ) = /o/2 and u±( 0 ) = 0 . Figure [5Jb) shows 
the W 2 dependence of 9 = arg(w^(0) • u + (T)) for several 
77 . Since the dynamical contribution is factored out, 9 
captures the geometrical contribution. For u >2 = 0, since 
Hk can be mapped to 9 = 0 for rj < 0.5, while 9 = n 
for rj > 0.5. On the other hand, for the finite 102 , which 
breaks the “chiral symmetry”, 6 deviates from 0 or n. In 
short, the quantization of the Hannay angle into 0 or n is 
protected by the “chiral symmetry”. Due to Eq. ([ 9 ]) , it is 
natural that the quantization condition for the Hannay 
angle is the same as for the Berry phase. 

The bulk-edge correspondence is useful in connecting 
the quantized Hannay angle and the topological character 
of the system in terms of edge states. In Fig. [3]jc), the 
frequency spectrum for the finite length system having 
the fixed boundary with LU 2 = 0 is plotted as a function 
of 77 . There is a clear signature of in-gap edge state dis¬ 
tinct from the bulk contribution for rj > 0.5 where 9 = ir, 
while there is no sign of edge state for 77 < 0.5 where 


9 = 0. Therefore, the quantized Hannay angle is used 
to detect the topological transition characterized by the 
appearance of localized edge modes. Here, we should be 
careful on the unit cell convention to establish the rela¬ 
tion between the Hannay angle and edge states, since the 
different unit cell conventions lead to the modified Hamil¬ 
tonian, which affects the Hannay angle m- In specific, 
the convention such that the boundary is in between the 
two neighboring unit cells should be employed. 

Figure [3^b) shows the frequency spectrum obtained 
with the free boundary condition. In contrast to the fixed 
boundary case, there is no sign of in-gap edge state irre¬ 
spective of the value of 77 . As discussed in the similar sit¬ 
uation [24], this contrast is actually the manifestation of 
the symmetry protection of the topological states. Recall 
that the quantization of the Hannay angle is protected 
by the “chiral symmetry”. For the fixed boundary con¬ 
dition, the “chiral symmetry” is preserved even after the 
boundary is introduced, since the diagonal elements of H 
are uniform even with the boundary. On the other hand, 
the free boundary breaks the “chiral symmetry”, since 
the diagonal elements of H becomes nonuniform due to 
the absence of a spring at the free boundary [23J • Then, 
the edge mode is not necessarily protected because of 
the symmetry breaking, even if the bulk parts share the 
same topological character for the two kinds of boundary 
conditions. 

To summarize, we derive a concise formula for the 
geometrical angle, i.e., the Hannay angle, in mechani¬ 
cal system by making use of the canonical transforma¬ 
tion that maps the Newtonian equation of motion to the 
Schrodinger type equation. Then, the use of the Hannay 
angle as a symmetry protected topological order param¬ 
eter is demonstrated using a simple spring-mass model. 
The importance of the symmetry in topological char¬ 
acterization is pointed out by analyzing the boundary 
condition dependence of the edge modes. It reveals yet 
another role of the symmetry in mechanical systems in 
relation to the bulk-edge correspondence. 
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Appendix: The Case with Hermitian f 


Let us start with the Hamiltonian 


H = p'iV'i + q'i^ i 3 q'r 


( 12 ) 
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where T is an Hermitian matrix. Now, we consider a 
canonical transformation using a generation function 


= -Wilijq'j + A+ iVZqaijq'j - mh 3 q 3 ■ ( 13 ) 


Here 7 = vf and 7 is hermitian under the condition that 
r is positive semidefinite. In the following, we assume 
this condition is fulfilled. Then, we have 


dF 

p ' i = M = + 


dF 


p\ = ^7 = -hijtfj + i Vhuqj, 
dF 

Pi = - = -i yfecfj'in + iqjiji, 

dF 

Pi = --qZ- = -iVZjijq’j +i7 ijQj, 


( 14 ) 

( 15 ) 

( 16 ) 
( 17 ) 


which gives us 


i Pi = ~ Qjlji), i Pi = ^(i Pi - 7 'ijQj), 

( 18 ) 

7 ij q j = V + ^). ^7ji = ^ + ipi )' 

( 19 ) 


and 


iPi = ^ (iPi + ^-7ji), i Pi = ^ (iPi + lij q'j) 1 

( 20 ) 


7 ij q j = ^('7*j^ - iPi), 9i7ji = ^(^7ji - iPi)- 


( 21 ) 


The Hamiltonian is transformed as 


H = - ^(}Pi- q ilil) (m - lij q j ) 

+ ^ (iP* + <7»7*z) (iP; + 7ijQj) 
=miijqj + miijPj- 


( 22 ) 


Then,the equations of motion are given as 

Qi = 17 ij q j, q i = njlju ( 23 ) 

Pi = -iPjlji, pi = -hijpj■ ( 24 ) 

If we choose an initial condition such that (<?■)* = < 7 ' 
and (jp'i)* = Pi, it is mapped to an initial condition such 
that (i Pi)* = 7 ^qj and (ip)* = qj'lji- Since the equation 
of motion implies that ipi and 7 i 3 q 3 , and i pi and q 3 T 3 i 
obey conjugate equations, respectively, such relationships 
are preserved in time evolution. On the other hand, the 
relations (ip,)* = 7 i3 q 3 and (ipi)* = q 3 ^ 3i imply (q()* = 
q[ and (p()* =p\. 
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